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ABSTRACT 

The results of the 1979 intensive field study to investigate continuous 
fumigation of the Nanticoke Generating Station's (NANTGS) plume are presented. 
Compared to a similar, but larger-scale study carried out in 1978 by the Atmospheric 
Environment Service, Ontario Hydro and the Ontario Ministry of the Environment, the 
fumigation zones were observed further out from the stacks and in the north-west 
quadrant. The intensity of fumigation (i.e., the ground level concentrations inside the 
fumigation zone) was generally lower in view of the relatively lower output of 
NANTGS during the study period. The data from the study are used to validate the 
Nanticoke Environmental Management Program (NEMP) fumigation model. The model 
is shown to agree with the observed data within the framework of the uncertainties in 
the various input parameters. 

Parallel studies to investigate the S0 2 oxidiation rate in the NANTGS 
plume, and the ambient levels of hydrocarbons near the TEXACO petroleum refinery 
are also discussed and recommendations are made for further study. The calculation 
of the SO^ oxidation rate was not possible because high sulphate levels found in 
samples collected upwind of the NANTGS made it difficult to correct the in-plume 
measurements for entrained sulphate. 

The ambient hydrocarbon monitoring shows no clear differences in the 
hydrocarbon concentrations measured upwind and down-wind of the refinery. 
However, the study provided an opportunity to investigate some new techniques in 
ambient hydrocarbon monitoring and in-plant measurements at a large refinery. 
Several shortcomings in the sampling methodology were recognized, and 
recommendations were made for future work in the Nanticoke area. 

Other aspects of the 1979 intensive field study are reported in detail in a 
number of companion reports — see References £-6 and 8. 
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1. Introduction 



The Air Resources Branch's Technology Development and Appraisal Section 
carried out an intensive field experimental program at Nanticoke during May 28 to 
June 14, 1979. The principal aim of this program was to characterise the dispersion of 
pollutants (particularly S0 2 ) from Ontario Hydro's Nanticoke Generating Station 
(NANTGS) under fumigation conditions, although parallel studies were also carried out 
to determine the SO ? oxidation rate in the plume, and also ambient hydrocarbon levels 
in the Nanticoke area (especially up - and downwind of Texaco), 

NANTGS is a coal fired electric power plant with a capacity of 4000 MW. 
Under present operating conditions, emissions are about 156,000 tonnes of SO ? per 
year. The gases are emitted through two 198 m stacks situated at the shoreline. The 
plumes from these stacks are generally very buoyant and under normal conditions are 
observed to rise to a height of about 400m from the stack base. 

Fumigation of NANTGS plumes occurs as a result of plume interaction with 
the thermal internal bounary layer (TIBL). The TIBL forms over land for gradient 
onshore/lake breeze flows. It is a convective boundary layer, capped by an inversion. 
The plumes are emitted into the stable layer and at some downwind location they 
intersect the TIBL causing fumigation, and relatively high ground level concentrations 
(gJc). 

The experiment was designed to obtain crosswind concentration profiles for 
each event, both inside the fumigation zone and close to the stack. This was 
supplemented by meteorological measurements for later interpretation of the data. 
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To measure the in-plume SO- concentrations an instrumented helicopter 
was used to collect airborne data, a ground-based mobile monitor to track the plume 
and to profile ground-level concentrations, and a mobile chemistry laboratory (the Air 
Monitoring Unit of the Monitoring and Instrumentation Unit of the Air Resources 
Branch) to measure half-hourly ground-level concentrations. Continuous 

meteorological measurements were obtained from the Ontario Hydro meteorological 
tower (temperature and winds at 10, 32 and 85M), and the Atmospheric Environment 
Serivces' (AES) weather station (hourly surface data). Minisonde and acoustic sounder 
data were used to define the thermal internal boundary layer (TIBL). The Atmospheric 
Environment Serivce also provided on-site forecasts throughout the study. 

While the helicopter was profiling the plume in the fumigation and near- 
stack areas, samples were also collected for obtaining the SO- oxidation rate in the 
plume, and measurements were made of ambient concentrations of several 
hydrocarbons, both upwind and downwind of the refinery. 

The rate of oxidation of SO- to SO. in the generating station plume under 
fumigation conditions is of special interest for two reasons: 

(1) under such conditions, a large volume of background air is mixed with 
the chimney emissions - a factor conducive to relatively high summertime oxidation 
rates (see, for example, reference 1); and 

(2) fumigations caused by plume interaction with a TIBL occur in the 
Nanticoke area with southerly flows. This means that the Hydro generating stations 
and Texaco refinery emissions can potentially be mixed, creating a highly oxidizing 
atmosphere of nitrogen oxides and hydrocarbons which could cause high S0 2 oxidation 
rates. 



The hydrocarbon sampling study was undertaken because the rate of SO- 
oxidation (and also phenomena such as ozone formation) depends on the type of 
hydrocarbons present in this mixture (olephins, for example, causing much higher 
oxidation rates than paraffins), as well as the concentrations. It should be pointed out 
that another objective for the hydrocarbon study was to test the usefulness of a 
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sampling method developed by the Monitoring and Instrumentation Development Unit 
of the Air Resources Branch, especially for estimating hydocarbon emissions downwind 
of the refinery (by computing the net flux across a vertical plane perpendicular to the 
mean wind direction). Although the Texaco refinery at Nanticoke was designed to 
keep hydorcarbon emissions to a minimum, an independent estimate of emission rates 
was considered to be necessary. At the time of the May/June 1979 field study, The 
Source Assessment Unit of the Air Resources Branch was making hydorcarbon 
concentration measurements in the headspace of. the various tanks in the Texaco tank 
farm, in order to estimate emission rates from this particular area of the refinery. 
Their experimental methods and results are described in the Air Resources Branch 
report TDA-25-8Q , but are also summarized in this report. 

This report presents the detailed data gathered by The Speical Studies Unit 
of the Air Resources Branch during the May/June 1979 study, and also summarizes the 
results obtained by the other participating groups (references for sources of detailed 
data obtained by these groups are given in the report). The interpretation of all the 
data, along with the validation of a mathematical model for plume fumigation in the 
Nanticoke area, is also included, together with recommendations for future work. 
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2, Experimental 

2.1 Meteorological Data 

2.1. a) Minisonde Data 

Minisondes were released hourly, and tracked using the double theodolite 
method, from two locations, one at the lakefront in the Peacock Point area, and the 
other 15-20 km, inland positioned according to wind direction. 

The azimuth and elevation angles of the balloon were measured every 30 
seconds by the theodolite operators, while the minisonde receiver provided a plot of 
temperature vs time. The temperature profiles were digitized by hand, and tabulated 
with the theodolite measurements for entry to a digital computer. 

2. Kb) Meteorological Forecasts 

General meteorological information was provided by The Atmospheric 
Environment Serivce, which operated an onsite meteorological forecast office for the 
duration of the study. In addition to providing weather forecasts throughout the day, 
many other parameters were determined including; gestrophic winds, temperature, 
dewpoint, solar radiation, wind speed and direction at 10 m and hourly observations of 
cloud cover. This data is summarized in AES's Interal Report SD-9-19 entitled 
"Summary of Synpotic Meteorological Conditions during Nanticoke Air Quality Field 
Study - May 28 - June 14, 1979". 

2.2 Ground Level Monitoring 

2.2. a) Ground Based Mobile Monitoring (SO-) 

A ground based mobile monitor was provided by Moniteq Ltd. The unit consisted 
of a vehicle with a COSPEC remote SO- sensor, TECO and SIGN-X ambient SO ? 
monitors. The unit was employed to find the plume, and to measure the ground level 
SO- concentrations and overhead burden in the fumigation area. The unit was in radio 
contact with the base station and relayed information regarding the position of the 
plume to the base, to be used for positioning the helicopter and Mobile Air Monitoring 
Unit (MAMU). 
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2.2,b) Ground Based Stationary Monitoring (SO-) 

The MAMU unit was used to measure ground level concentrations of SO- 
along with a number of other parameters (including meteorological data) at a 
stationary point. The MAMU unit was also in communication with, and positioned by, 
the base station. This data, including half hourly S0 2 concentrations at various 
locations, may be found in TDA report #TDA 03-80. 

2.2.c) Ground Level Hydrocarbon Monitoring 

To determine the background hydrocarbon concentrations and also any 

contributions from the refinery and other local sources, a preliminary study was 

3 

carried out by the Ontario Research Foundation to develop a methodology to monitor 

ambient levels of hydrocarbons. The methodology consisted of 1) collecting ambient 
air samples in five-layer (Mylar, aluminum, Mylar, teflon, and polyester) bags obtained 
from Calibrated Instruments Inc.; 2) desorbing the concentrated sample into a gas 
chromatograph equipped with a flame ionization detector. This method, however, 
met with little success due to problems encountered in preserving ambient-level 
samples in bags, and design problems with the desorption unit. Subsequently, a method 
to determine the concentrations of various hydrocarbons using direct adsorption onto a 
pre-concentrating cartridge was developed by the Monitoring and Instrumentation Unit 
of the Air Resources Branch, 

The collection and anlaysis of the hydrocarbon samples was primarly the 
responsibility of the Monitoring and Instrumentation Development Unit, and the 
experimental techniques and complete analysis results are described in TDA Report 
No. 03-80 . Briefly, the ambient air was sampled through a glass cartridge packed 
with a Molecular Sieve 5A adsorbing medium for 20 to 60 minutes at a constant, 
measured flow rate. For the samples collected using the direct injection method, the 
sample was drawn over the adsorbent from the mobile monitoring vehicle's sampling 
manifold for 30 minutes, after which it was automatically injected into the gas 
chromatograph (H.P. 5830A gas chromatograph with a 9 ft., 1/8" O.D. stainless steel 
column packed with Porapak N 80 - 100 mesh). The components analyzed for are 
indicated in Figures 1 - W. 
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The ground based portable samples were collected by drawing air through a 
cartridge for 60 minutes at a flow rate of 50 ml/min, Nutech samplers were used, 
with restricting orifices limiting the flow to 50 ml/min. The ground level cartridges 
were located up- or down-wind of the plant at least 100 feet from the portable 
generator to avoid sampling the generator exhust. Handing and setup blanks were also 
analyzed, and all levels reported here were corrected for blank values. 

2.3 Airborne Monitoring 
2.3.a) S0 2 Profiling 

The airborne SO- concentrations were measured using a Sign-X SO ? 
analyser with a sampling nozzle mounted on the landing skid, extending in front of the 
helicopter to minimize the effects of the downwash from the chopper blades. The 
Sign-X was calibrated before, and several times during the field study. It was found 
to have a linear response, and on days that calibration did not take place the average 
of the slopes of the calibration curves was used. The slopes of the calibration curves 
varied by less than 10%. 

Each flight consisted of a set of crosswind passes in the fumigation area folio wed 
by a set near the stack. A background flight was carried out on most days by sampling 
at plume height well upwind of the stack. The fumigation area was found by visual 
observation followed by a low pass in the helicopter. Upon finding the fumigation 
area, a flight path was established perpendicular to the plume axis and sampling 
began. The first pass was made at the low altitude followed by the middle, then the 
upper in an »S' shaped path to avoid interference from the helicopter. A typical set of 
passes is shown below: 

900 ft 

600 ft Repeat 

300 ft 

As many passes as possible were completed at one downwind distance. However, 
the fumigation area often shifted, making it necessary to move further downwind to 
get back into the fumigation area, thus cutting down on the number of passes that 
were used in forming the crosswind concentration average (see below). 
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The pilot was instructed to maintain a straight path and a constant 
airspeed of 30 mph. This, however often varied from 20-40 mph and the bearing 
sometimes changed by as much as 30 degrees during a given pass. The plume 
near the stack was sampled using the same '$' shaped pattern, but at altitudes 
dictated by the plume height and vertical dispersion. 

The helicopter SO^ profile data were positioned such that all travreses in 
a set started from a common point and went in the same direction. The 
traverses and starting points are marked on the map, and summarized in a table 
for each day. The data are grouped so that all taverses at one downwind 
distance are contained in a "run 1 . Each flight consists of 2-4 runs at different 
sampling locations. 

2,3.b) SO- Oxidation 

The data for the determination of the oxidation rate of SO,* were 
collected using 47 mm filter packs mounted in the helicopter. The sampling 
probe consisted of ah isokinetic teflon nozzle mounted on a landing skid 
extending in front of the helicopter to minimize the effects of the downwash 
from the blades. The nozzle was connected with ft" I.D. smooth walled TEFLON 
tubing to a two stage polyethylene SWINNEX filter pack, The filter pack 
consisted of two filters in series, the first a Whatman 40 prefilter to trap 
particulates, followed by a Whatman 41 filter impregnated with potassium 
carbonate-glycerol solution to trap gaseous SO^. The air flow rates through the 
filter packs were monitored with a calibrated BROOKS rotameter and found to 
be quite constant ( less than 10% variation). 

The filters were exposed during the flights that were designed to 
provide SO~ profiling data. Samples were collected at two locations downwind 
of the stack, as well as one upwind from the stack (to collect a background 
sample). Section 2.3a describes the flight path, and the individual flights are 
described in Appendix II. A total of six to fifteen traverses was completed in the 
fumigation area over 20-60 minutes, with the pumps left on during turns (outside 
the plume). Upon completion of the flights in the fumigation area, the filter pack 
was changed (and stored in an air tight container) and the near-stack sample was 
collected in a similar fashion at altitudes dictated by plume height and vertical 
spread, usually 1000-2000 feet above ground level. The background sample was 
collected by making several traverses (over 30 min) well upwind of the 
generating station at plume height. 
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The impregnation, and chemical analysis of the filters was carried 
out by Ontario Hydro. The Whatman 40 prefilters were extracted using two 15 
ml portions of hot deionized water, and the Whatman II (impregnated) filters 
were extracted using two 25 ml portions of 0.03% hydrogen peroxide solution. 
The extracted solutions were then analyzed using ion chromatography. The 
results in each case were reported as micrograms of sulphate per filter. 

2,3.c) Airborne Hydrocarbon Sampling 

The airborne hydrocarbon samples were collected by helicopter over 
various time periods while traversing up- or down-wind of the petroleum 
refinery. The sample was collected through a V stainless steel sample line at a 
flow rate of about 50 ml rnin (monitored with a Hastings Mass Flowmeter) 
through the cartridge. Care was taken to avoid sampling helicopter exhaust 
fumes, by turning the helicopter into the wind and shutting the pumps off during 
turns. Test flights were performed specifically to catch and obtain a 
'fingerprint* of the helicopter exhaust (see below). The cartridges were then 
thermally desorbed into a gas chromatograph for analysis. 

2.4. Source Assessment Measurement of the Tank Farm 

The Source Assessment Unit of the Air Resources Branch was also 
making concentration measurements In the headspace of the various tanks at the 
tank farm (which were suspected to be a major source of emissions) during the 

time the ambient sampling was taking place. The experimental methods and 

5 
results are described in Air Resources Branch/TDA report #25-80 . Briefly, 

evacuated flasks were used to collect samples of the headspace gas, which were 

then delivered to the laboratory for G.C. analysis. 
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3. Results and Discussion 
3.1 Meteorological Data 

The minisonde data presented in Appendix IN consists of pairs of 
temperature profiles for each balloon flight (usually 1 per hour on the hour) 
throughout each day. The M5-1 profile is from the lakefront site (at Peacock 
Point), while the MS-2 profile is from the inland site (see map in Appendix II for 
the location of this site, which varied from day to day). The shape of the 
thermal internal boundary layer (TIBL) is assumed to be parabolic (see Appendix 
III) and given by: 

L(x) = AX 5 * 

where A is a constant, L is the height of the boundary layer at a downwind 
distance x. The value of A for input to the model is calculated by knowing L at a 
given downwind distance, which is defined as the distance from the shoreline to 
the point of measurement in a direction parallel to the mean wind. These data 
are used as input parameters to the NEMP fumigation model which is discussed 
in Appendix I. 

3.2 Ground Level SO- Monitoring 

The ground level SO monitoring was carried out by two vehicles; the 
Monitoring and Instrumentation Unit's Mobile Air Monitoring Unit (MAMU) for 
stationary measurements of ground level SO , and a traversing vehicle to track 
the plume with a remote SO- analyser and measure ground level SO- in the 
fumigation zone. 

The results of the MAMU survey are presented in ARB-TDA Report 
03-80 . Some of the half hourly average SO- concentrations have been used in 
this study for comparison to the ground level SO- concentrations predicted by 
the NEMP fumigation model (Appendix I), The half hourly SO- average 
concentrations have been extracted and tabulated in Table 1, While examining 
the data it is important to note that data presented are half-hourly averages (not 
hourly averages as required for comparison to the Provincial Ambient Air 
Quality criterion) collected during a period that the NANTGS was operating at a 
reduced rate (see below). 
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The data from the ground based profiling may be found in Moniteq 
report MTR80-15 , The profiles of ground level SCU concentration in the 
fumigation zone and overhead burden of SO« are presented as functions of cross- 
wind distance. Information on sampling times and locations during this project is 
summarized in Table 2. 



The ground level SO- profiles were plotted and averaged to give an 
averaged profile of ground level SG ? concentrations as a function of crosswind 
distance in the fumigation zone. These averaged profiles were compared to 
similar averaged profiles for helicopter data collected at 300 ft. to ensure that 
the model's assumption of uniform vertical mixing within the TIBL is valid. In 
general the ground level concentrations are within 10% of the concentrations at 
300 feet above ground (See Appendix I, Table 3). 

The averaged ground level SO- prolfiles were also used for comparison to 
predicted levels calculated by the NEMP fumigation model (see Appendix I for a 
detailed discussion), 

3.3 Airborne Monitoring 

3.3.a Airborne SCU Monitoring 

Table 3 contains information showing the flight/run indentification, the 
time and location of sampling and the number of crosswind transects completed 
at each downwind distance. The actual data for each day are contained in 
Appendix II. Included for each day is a description and evaluation of the day's 
sampling, a table summarizing the locations and starting points for each group of 
transects, a map showing the location of the helicopter sampling, the MAMU unit 
and the minsonde release sites, and the individual plots of SO- concentration as a 
function of crosswind distance at each altitude. 

A total of 9 sets of traverses were completed within the fumigation zone 
during the 8 sampling days of the study period. The fumigation zone was 
sampled at various distances from the NANTGS ranging from 12-30 km from the 
source under various operating conditions. The NANTGS net megawatt output 
and estimated SO- emissions are summarized in Table 4. 

The SO- airborne data has been used to validate the NEMP fumigation 
model and to investigate the fumigation of the NANTGS plume. The validation 
of the model is described in detail in Appendix I. 
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Briefly, the NEMP model predicts the ground level concentrations of 
pollutants emitted from a tall stack into a shoreline environment under 
fumigation conditions. The model assumes that the dispersion of the plume 
proceeds in two steps: 1) The plume is assumed to disperse in the standard 
Gaussian manner in the stable layer above the TIBL, followed by 2) dispersion 
into the TIBL. The dispersion into the TIBL is treated as an elevated area source 
located at the under-surface of the TIBL envelope. 

The near-stack elevated SO- profiles have been used to verify the 
assumption of standard Gaussian dispersion in the stable layer, while the SO- 
profiles in the fumigation zone have been used along with various ground level 
data for comparison to ground level concentrations predicted by the model. In 
general, the point by point comparisons of the experimental data to the 
predicted cross-wind profiles agree very well near the centre of mass, but the 
predicted profiles tend to fall off more rapidly beyond 500 m from the centre of 
mass than the experimental data. However, the relatively low concentrations 
measured in the tail ends may have a larger margin of error thus making the 
comparison of predicted concentrations to measured values difficult in the tails 
ends of the profiles. The scattergram of the observed to predicted levels are 
seen to be evenly distributed about the one-to-one line (Appendix I, Fig. 2). 

3.3.b Airborne Measurement of SO- Oxidation Rate 

I 

Since the data for this study of the SO- — ^ SO. oxidation rate were 
collected as part of a plume fumigation study, the general meteorolgical 
conditions consisted of on-shore flows throughout the sampling period. The 
weight of particulate sulphate trapped on the impregnated filters ranged from 
120-1000 ug per filter (as sulphate). Blanks from each batch of filters were also 
analyzed, and blank values were subtracted from the reported values. There 
were some obvious cases of contamination of blank filters, and in these cases an 
average blank value was subtracted (2.0 ug/filter for sulphate, and 1^.0 ug/filter 
for SO-). Handling and setup blanks were also analysed, and values were 
compared to true blanks with no significant discrepancies observed. 

The data including location and time of sampling are summarized in Table 
5. The amounts of SO- and sulphate are reported as micrograms of sulphur per 
cubic meter of air sampled. A daily flight description is presented in Appendix 
II. 

To calculate the SO- oxidation rates it is necessary to correct the SO- and 
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total sulphate values in the plume samples for the contributions from entrained 
ambient air. However, upon examination of the background sulphate values, it 
can be seen that although the SO^ values are reasonable, the sulphate levels in 
the background sample are often greater than those found in the sample 
collected in the fumigation area, which makes such a correction impossible. 
These background samples were collected at 1000 ft. a.g.l., while the fumigation 
samples were collected from 300-900 ft. a.g.l. The background sulphate 
measurements were compared to ground based Hi-volume sampler data collected 
by Ontario Hydro (Peacock Point) and MOE (5 samplers in the Nanticoke region) 
on the same days, and were found to be generally higher (see also Hydro results 
- Table 5). This suggests that there may have been a stratified distribution of 
background sulphate with altitude on the sampling days, possibly due to emissions 
from elevated sources on the U.S. side of Lake Erie, which were decoupled from 
the surface while being carried over the lake. 

Further evidence for the possiblity of such a stratification may be found in 

the data collected during the Atmospheric Environment Service Nanticoke 

2 
Shoreline Diffusion Experiment of 1978 (Anlauf et al) . Aircarft measurements 

made on June 7 (a day with strong southerly flows) show a very high light 

scattering coefficient, increasing with altitude through the mixing layer. The 

profile of the light scattering measurements is muitilayered, suggesting the 

possiblity of stratification of sulphate with height (since sulfate is strongly 

correlated with the particles in the light scattering size range). If such a 

stratification exists it is possible that the ambient air entrained into the plume 

in the fumigation area (sampled at altitudes less than 1000 ft.) has less SO. than 

the "background" air collected above 1000 ft., upwind of the stacks, and the 

sulphate concentration in the so-called "background" sample may be higher than 

that in the mixture of plume and entrained air in the fumigation zone. Thus the 

present results make it very difficult to correct the in-plume measurements for 

entrained sulphate, and the SO ? oxidation rate cannot be calculated. 

3.4 Hydrocarbon measurements - Airborne and Ground Level 

When analysing and interpreting the hydrocarbon data, it must be 
remembered that these studies were of a preliminary nature, and sampling and 
analysis methods were not optimal (see below). Therefore, the conclusions 
should be regarded as tentative only. In future work, an attempt will be made to 
improve on the measurement methodology. 
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For each of the sampling days the downwind measurements have been 
averaged, and compared to the average of the upwind values. These average 
concentrations are plotted on Figures 1 to 10, one for each day's sampling and a 
map showing the sampling locations is given on pp. 42-50. The number of values 
used in computing each average is shown in brackets in the legend for these 
Figures. On some of the days, upwind samples were not collected, and only 
downwind values appear. Airborne samples were collected at various altitudes, 
up- and down-wind on two days during the study (Figures 1 to 9). On both days, 
high levels of ethene and 1,3-butadiene were detected in both up-and down-wind 
samples. Both of these are products of incomplete combustion, and were found 
in the helicopter exhaust (see below). This suggests that in spite of the 
precautions taken to avoid sampling exhaust fumes, some contamination may 
have taken place, although some of the other components found in the exhaust 
(n-hexane, n-heptane) did not appear. 

Many of the ground based samples (field cartridges and direct injection 
method) also show high levels of ethene and 1,3-butadiene in both the up- and 
down-wind samples. The levels of the other components show no upwind - 
downwind differences when averaged. However, there are some individual down- 
wind cartridges exhibiting very high levels of some of the higher molecular 
weight components (see Figure 5 and Figure 8). These high concentration 
samples may have sampled "puffs" of hydrocarbons that were emitted over a 
short period of time, and were not picked up on other cartridges. On other days 
reasonably high levels of n-pentane, 3 -methylbutene and n-hexane were 
detected down-wind of the refinery (on Concession Road No. 3, 2 km east of 
Nanticoke Road), but no up-wind samples were collected for comparison (Figures 
6 and 7). 

Several cartridges were exposed to obtain 'fingerprints' of various sources 
including helicopter exhaust, abandoned gas wells in the area, pipelines and 
underwater gas wells at the lakeshore. These fingerprints are illustrated in 
Figures 11 to 14. The helicopter exhaust can be seen to contain almost all the 
hydrocarbons found in ambient air samples, and especially ethene and 1,3- 
butadiene. The sample concentrations collected near leaking gas wells show high 
levels of i-butane, n-butane/butane and propane/propene, and lower levels of 
some of the heavier components. The sample concentrations collected at Long 
Point to determine the hydrocarbons in air coming from Lake Erie were low level 
with the exception of 1,3-butadiene which was found at 30 ppb. The cartridges 
exposed downwind of the area in which the Texaco pipelines cross Concession 
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Road No,2 were found to have quite high levels of some of the heavier 
hydrocarbons (particularly toluene, benzene and hexane). However, after 
discussions with Texaco personnel (R. Cameron, March 3, 1980) it was found 
that these samples may have been collected near a sewer leading to the 
refinery's water treatment plants, on which a man hole had a cover that was 
improperly sealed at the time of sampling. 

The measurements made by the Source Assessment Unit at the tank farm 
provided a 'fingerprint' for the components in the headspace of the various tanks, 
but these could not be related to the observed ambient measurements. However, 
they did provide data regarding the constituents in the headspace of the various 
tanks examined. In particular, it is interesting to note that several 
photochemically reactive species, especially propylene and i-butylene, were 
found to occur at significant levels in the headspace of some of the tanks (see 
Table 6). 

In general, there were a number of problems associated with both the 
ambient and in-plant sampling and analysis methods. 

In the ambient sampling methodology some of the problems encountered 
were: 

1) difficulty in locating the hydrocarbon plume due to the variability of the 
winds around the refinery and also, the probable low hydrocarbon emission rates. 
The only wind data available was from the Ontario Hydro meteorological station 
located at Jarvis, and MAMU mobile laboratory which was often far from the 
site of the measurement; 

2) possible sources of sample contamination due to the products of 
combustion, particularly with respect to the airborne sampling; 

3) possible loss of samples from cartridges prior to G.C. analysis during 
transfer to the thermal desorber; 

4) the possibility of ambiguous identification of components due to similar 
retention times on the G.C. (E. Singer, personal communication). 

The Source Assessment Unit also experienced some difficulties in their 
experimental method including: 

1) loss of sample due to condensation problems in the evacuated flask; 

2) losses in the unheated line running from tank to sampler; 

3) possible ambiguous analysis of some compounds due to similar G.C. 

retention times. 
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(f. Conclusions and Recommendations 



4.1 Conclusions 



The characteristics of dispersion of NANTGS plumes under fumigation 
conditions are presented. Comparing with the results of previous studies it is 
clear that the TIBL at Nanticoke is highly variable. Consequently the downwind 
location of the fumigation zone is also highly variable. It is apparent, however, 
that many fumigation episodes at Nanticoke may be occurring beyond a 
downwind distance of 15 km. This is attributed to the formation of a rather 
shallow TIBL at Nanticoke 



The short range fumigation model to predict downwind glc has been shown 
to agree with the experimental data. Similar agreement was also observed for 
the data set of the 1978 AES study. This shows that the model is able to 
reproduce the physical system of continuous fumigation well. 

The calculation of the SO- oxidation rate from these data was not possible 
because the in-plume measurements could not be corrected for entrained 
sulphate. 

This study has provided an opportunity to field test some new techniques in 
ambient hydrocarbon monitoring and in-plant measurements of hydrocarbons at a 
petroleum refinery. Further tests of the sampling method and field techniques 
will be required to perfect the technique. However, the ambient sampling 
results show no clear differences in individual hydrocarbon concentrations up- 
wind and down-wind of the refinery. These observations may be compared to 

7 
some measurements reported by Westberg et al ., where a definite refinery 

effect on ambient hydrocarbon concentrations could be detected as far as 25 km 
down-wind. The in-plant measurements, in spite of the difficulties encountered 
with the sampling methods, provided some useful qualitative data on the 
constituents of the headspace gases in the various tanks in the refinery. 



- 17- 



k.2 Recommendations 

Further study is required to investigate the occurence of fumigation of the 
NANTGS plume at downwind distances greater than 15 km. 

Future studies of the SO. oxidation rate in Nanticoke should include an 
investigation into the possible stratification of sulphate concentration with 
height in the ambient air, particularly when there are southerly flows. To make 
meaningful background corrections, a background sample must be representative 
of the ambient air entrained in the plume, and collection should be carried out at 
different altitudes. The inclusion of a nephelometer in the instrument package, 
as a real-time surogate for sulphate concentrataion measurements, is strongly 
recommended for any future studies of this type. 

The following improvements are recommended for futher studies of 
ambient hydrocarbon levels near the Texaco refinery. 

1) Since vertical concentration measurements also are required to properly 
determine a hydrocarbon flux profile and sample contamination is difficult to 
avoid using a helicopter, alternative airborne sampling platforms should be 
investigated. 

2) Meteorological measurements should be made at the sampling site, to 
assist in accurately locating the samplers within the plume. 

3) An attempt should be made to assist in the location of the sampling 
cartridges by real-time measurements of the hydrocarbon plume from the 
refinery (possibly by using a total hydrocarbon analyser mounted in a traversing 
vehicle). 

%) Duplicate ambient samples should be collected to check the precision of 
the method. 

5) Further design and development work should be undertaken to improve 
sampling and analysis methods for both the ambient and the in-plant 
measurements. 

These recommendations will be incorporated into a second hydrocarbon 
study slated to take place in summer or fall of 1980. 
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Table 1 
Average Half Hourly SO,* Concentrations from MAMU 



Date 


Time 


Distance Maximum Average SO ? 
From NANTGS 30 min average concentration 










(km) SO„ concentration (ppm) for Monitoring 






2 


period (ppm) 


May 28 


1207-1355 


2 


.009 


003 


29 


1048-1254 


16 


.014 


.005 


29 


1341-1501 


11 


.052 


025 


30 


1036-1351 


10.5 


.015 


,011 


30 


1427-1515 


U 


.050 


043 


30 


1533-1630 


16 


.024 


017 


30 


1650-1741 


18 


.035 


027 


30 


1410-1928 


10.5 


.035 


009 


31 


1612-1927 


10.5 


.012 


003 


31 


1931-1406 


10.5 


.059 


013 


3un. 2 


1710-1813 


15 


.091 


057 


3 


1628-1731 


5.5 


.066 


045 


% 


1133-1327 


13.5 


.016 


013 


ft 


1350-1459 


15.5 


.130 


081 


4 


1648-1742 


16.5 


.140 


101 


6 


1539-1642 


15 


0$0 


068> 


6 


17^9-1044 


10.5 


056 


018 


7 


1217-1317 


16 


067 


046 


7 


1358-1504 


18 


190 


123 


7 


1613-1749 


27 


160 


123 


S 


0531-0628 


10.5 


170 


154 


8 


1341-1456 


17 


031 


022 


S 


2138-0828 


10.5 


002 


001 


9 


1037-1234 


1.5 


019 


.012 


9 


1331-906 


10.5 


014 


007 


10 


1014-0829 


10.5 


Oil 


002 


11 


1129-0554 


10.5 


003 


002 


12 


1011-1126 


2 


001 


001 


12 


1402-1517 


1.5 


001 


001 


12 


1534-1716 


1.5 


.001 


001 


13 


1458-1558 


16 


090 


.077 


13 


1642-1818 


23 


140 


.113 


14 


1438-1522 


16 


330 


244 


14 


1714-1801 


23 


180 


157 
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Table 2 
Summary of Ground Level Plume Profiling 



DATE 



TIME 



LOCATION DISTANCE NUMBER OF 

from NANTGS CROSSWIND 
(KM) TRAVESES 



May 30 



1505-1545 
1550-1735 



Cone. 12 
Hwy,40 



17 
16 



May 31 



1112-1140 Hwy. 6, Port Dover 



10.5 



Jun. 4 



1300-1435 


Hwy. 40 


16 


6 


1525-1600 


Hwy, 3 


18 


2 


1600-1632 


Hwy. 40 


16 


3 


1640-1720 


Rainham Rd. 


18 


5 



Jun. 6 



1 500- 1611 Rd. 5- Port Dover 



14.3 



Jun, 7 



1051-1142 



Hwy. 40 
Hwy, 54 



16 
30 



Jun. 14 



1242-1525 
1608-1728 



Hwy 40 
Hwy. 54 



16 
30 



21 
8 
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TABLE 3 
Summary of Helicopter SO- Profiling 



DATE 


TIME 


FLIGHT/RUN 


LOCATION 


DISTANCE FROM 

NANTGS 

KM 


NUMBER OF 
S0 2 PROFILES 


MAY 29 
May 29 


1443-1525 
1537-1608 


5/1 
5/2 


1 km SW of HWY 40 
Cone. 1 


ABORTED 




May 30 
May 30 


1735-1525 
1828-1845 


6/1 
6/2 


HWY.40 
Cone. 1 


16.5 

2.0 


6 
6 


May 31 


1328-1415 
1432-1452 


7/1 
7/2 


Turkey Point 
Stelco Pier 


30 
3.7 


6 
9 


Jun. 4 
Jun. 4 


1357-1430 
1458-1505 


8/1 
8/3 


Nelles Corner 
Cone. 1 


19 
3.0 


5 
3 


Jun. 4 
Jun. 4 
Jun. 4 


1635-1656 
1709-1719 
1743-1751 


9/1 
9/2 
9/3 


Grand River 

East River 

Cone. 1 


2.6 
3.2 

3.0 


D 

3 
4 


Jun. 6 
Jun. 6 


1445-1543 
1555-1620 


10/1 
10/2 


Port Dover 
Stelco 


12.3 
3.0 


10 
9 


Jun. 7 
Jun. 7 
Jun. 7 


1315-1332 
1350-1420 
1432-1455 


12/1 
12/2 
12/3 


Sandusk Rd. 

HWY. 40 

10 km NE of HWY 40 


2.4 

16.7 

25 


7 
7 
6 


Jun. 13 
Jun, 13 
Jun. 13 


1605-1641 
1652-1716 


18/1 
18/2 
18/3 


Aborted 

HWY 40 

Sandusk Rd. 


— 
17.8 
3.6 


9 
11 


Jun. 14 
Jun. 14 
Jun. 14 


1528-1350 
1400-1418 
1428-1447 


20/1 
20/2 
20/3 


HWY 40 

Sandusk Rd. 

upwind 


16.7 
3.6 


15 

9 


Jun. 14 
Jun. 14 


1545-1640 
1655-1709 


21/1 
21/2 


5 km NE HWY 40 
Sandusk Rd. 


21.3 
3.6 


12 
9 
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Table 4 
NANTGS Hourly Emission Rates 



Date Hour Total Output Estimated SO Emitted* 

Mgrf 1 

May 29 11 475 7.3 

7.3 
7.3 
7.3 
7.3 
7.3 
7.3 
7.3 

May 30 11 471 7.2 

7.2 
7.2 
7.2 
7.2 
7.3 
7.4 
7.4 
3.8 

May 31 11 482 8.0 

8.0 
8.0 
8.0 
8.0 
7.9 
7.9 
7.4 

22.8 
22.7 
27.6 

22.5 
22.5 
22.5 
20.9 



four 


Total Output 




MW 


11 


475 


12 


476 


13 


476 


14 


476 


15 


476 


16 


477 


17 


476 


18 


477 


11 


471 


12 


472 


13 


472 


14 


472 


15 


469 


16 


473 


17 


478 


18 


428 


19 


247 


11 


482 


12 


483 


13 


485 


14 


483 


15 


486 


16 


481 


17 


480 


18 


449 


11 


1441 


12 


1437 


13 


1433 


14 


1426 


15 


1423 


16 


1422 


17 


1424 


18 


1323 
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Table I (Cont'd) 

Date Hour Total Output Estimated SO- Emitted 

MW Mg hi 



Jun. 6 11 986 15.3 

15.3 
15.3 
15.2 
15.2 
15.2 
15.2 

18.2 
16.3 

IM 
16.4 
16.4 
16.4 
16.4 
16.4 
16.4 

hii i ii 16.6 

16.7 
16.7 
16.8 
15.1 
15.1 
12.9 
^.5 
3.7 



Jun. 14 11 963 15.9 

15.9 
15.9 
15.9 
15.9 
15.8 
13.0 
7.5 



Calculated by: MW, tota i\ .h x multiplier 

where multiplier = weekly gross/net MW ratio x 156 Mg coal/h x %S x 2(S to SO. ) 

500 MWh *" 



(See letter from D.K. Gillies, Ontario Hydro, to G. Van Volkenburgh, MOE, dated 
June 5, 1980, File NE-OS-07-02.) 
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11 


986 


12 


966 


13 


969 


14 


965 


15 


965 


16 


961 


17 


965 


11 


1064 


12 


953 


13 


959 


14 


963 


15 


963 


16 


962 


17 


959 


18 


960 


19 


960 


11 


955 


12 


958 


13 


960 


14 


963 


15 


869 


16 


867 


17 


743 


18 


486 


19 


213 


11 


963 


12 


963 


13 


963 


14 


964 


15 


962 


17 


960 


18 


786 


19 


452 



Table 4b 
Sulphur Content in Coal 



Week of 


Coal Analysis 




Sulphur Content as Fired 




% S 


May 30/79 


2.24 


June 6 


2.41 


June 13 


2.61 


June 20 


2.47 
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Table 5 



Summary of Oxidation Rate Filter Analysis Results 



May 30 



May 31 



June 4 






June 6 



June 7 



June 13 



June 14 



,T/Run 


Tim< 


* of 


Distance 


Volume 


SG^ 


so 2 


OH PEACOCK 


NEMP SO u 




Sampling 


From Stack 


Samaled 
(m T ) 


PT SO. 
Result! 
ug S/m 


RESULTS 








(km) 


ug S/m 


ug S/m 


(Ave. of 
















5 Stations) 


5/1 


1445- 


■1530 


19,0 


1.01 


2,14 


51.8 






5/2 


1537- 


■1611 


2.3 


.75 


3,0 


331 






6/1 


1735- 


■1815 


16,6 


.94 


3.4 


53.5 






6/2 


1828- 


•1848 


2,0 


,46 


1.09 


146 




2.36 


7/1 


1328-1409 


30.1 


1.0 


1.67 


52.3 






7/2 


1432- 


•1454 


3.7 


,53 


5.72 


358 






7/3 


1525- 


•1555 


Background 


.72 


3.00 


19.8 


2.3 




S/l 


1357- 


•1435 


18,8 


.93 


2.06 


119 


- 




8/3 


1458- 


•1509 


2.9 


.30 


2.78 


559 






9/1 


1635- 


■1600 


25.7 


.58 


3.3 


119 






9/2 


1709- 


-1724 


32.8 


.37 


3,15 


128 






9/3 


1743-1754 


2.9 


.26 


1.28 


344 






10/1 


1445- 


-1550 


12.3 


1.60 


3.33 


87 






10/2 


1555- 


-1623 


2.7 


.65 


9.23 


806 






11/1 


1725- 


■ 1755 


Background 


.70 


5.38 


13.3 






12/1 


1315- 


•1334 


2.4 


.44 


6.81 


203 






12/2 


1350- 


-1425 


16.7 


.80 


6.87 


112 






12/3 


1432- 


-1500 


25.1 


.64 


7.97 


160 






13/1 


1600- 


-1630 


Background 


.71 


6.60 


^.5 


9.2 




18/2 


1605- 


-1645 


17,8 


.99 


.33 


51.8 






18/3 


1652- 


■1720 


3.6 


.70 


- 


287 






19/1 


1804- 


-1830 


Background 


.63 


1.32 


9.0 






20/1 


1258- 


-1355 


16,7 


1.4 


1,66 








20/2 


1400- 


-1420 


3.6 


.50 


2.0 


282 






20/3 


1424- 


-1447 


Background 


.60 


3.33 


14.5 


1.5 




21/1 


1545- 


-1645 


21.3 


1.50 


3.0 


65.1 




2.2 


21/2 


1655-1711 


2,6 


.40 


5.0 


262 




2.2 



Table 6 

Analysis of Headspace Gases 

(From TDA 25-80 A Preliminary Study of the 

Emission of Hydrocarbons from TEXACO Nanticoke Refinery) 



Tank 
Content 


Methane 
(ppm) 


Ethyl- 
ene 
(ppm) 


Ethane 
(ppm) 


Propane/ 

Propylene 

(ppm) 


Iso- 

butylene 

(ppm) 


N-butane/ 

butadiene 

(ppm) 


Benzene 
(ppm) 


Toluene 
(ppm) 


Ethyl- 
benzene 
(ppm) 


Xylene 
(ppm) 


Styrei 
(ppm) 


crude 


1.2 


nd 


6.6 


12.6 


4.2 


nd 


nd 


nd 


nd 


nd 


nd 


gasoline 


0.7 


0.2 


0.3 


3.0 


389.7 


94.2 


8.4 


5.2 


0.41 


12.7 


0.9 


jet fuel 


us 


0.04 


0.1 


6.7 


19.1 


nd 


- 


29.3 


14.7 


33.7 


2.6 


bunker 


409 


40 


661 


323 


189 


95 


_ 


68 


15 


62 
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map 1 

May 28, 1979 
Hydrocarbon Sampling Locations 
A Ground level (downwind) 
B Ground level (upwind) 
C Airborne (dowro-riid) 
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Hydrocarbon Sanpling Location 
B Ground level upwind 
C Ground level downwind (2) 
E Ground level downwind (2) 
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MAP 3 

May 30, 1979 
Hydrocarbon Sampling Locations 

A Ground level upwind ,-,) 

B Ground level downwind (1) 
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■ map 4 
June 3, 1979 
Hydrocarbon Sarrpling Locations 

A Ground level upwind (1) 
B Ground level downwind (3) 
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Hydrocarbon Sampling Locations 

A Ground level upwind (1) 
B Ground level downwind (3) 
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MAP 6 
June 8, 1979 
Hydrocarbon Sampling Locations 
A Ground level downwind (2) 
B Ground level downwind (2) 
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MAP 7 
June 10, 1979 
Hydrocarbon Sampling Locations 

A Ground level downwind (7) 
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MAP 8 

June 12, 1979 

Hydrocarbon Sampling Locations 

G Ground level downwind (7) 

H Ground level downwind (4) 

K Ground level downwind (5) 

S Airborne Upwind (2) 

T Airborne Downwind (3) 
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Locations 

A Ground level upwind (7) 
B Ground level downwind (5) 
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Appendix I - N.E.M.P. Fumigation Model Validation 



i. DISCUSSION OF THE MODEL : 

Dispersion in the stable layer is assumed in the model to proceed in the 
standard gaussian manner. Since the downwind distance between the source and 
the fumigation zone on most occasions is relatively short effects of wind shear on 
the crosswind distribution of the plume in the stable layer are small. Therefore, a 
Gaussian crosswind distribution of the plume in this layer is not a very restrictive 
assumption. A Gaussian vertical distribution, however, may be unrealistic. 
However, as shown later, from the form of the final equation used to compute the 
glc the exact shape of the plume in the vertical direction in the stable layer is not 
very critical to the derived results. 

Dispersion into the TIBL is assumed to proceed from an elevated area source 
co-located with the under surface of the TIBL envelope (height surface). The 
strength of this area source is derived from the entrainment of the stable-layer 
plume by the rate of rise of TIBL height and injection of this plume into the TIBL 
by processes which disperse it in the stable layer. It is assumed that pollutants 
released into the TIBL at its top have lost their buoyancy and disperse quickly into 
it as passive contaminants. The crosswind spread with downwind distance inside 
the TIBL of particles released at its top are, therefore, prescribed in a manner 
similar to Lamb's (1979). This leads to the following equation for the glc within the 
TIBL: 



C6*r/) \JSl 



where 
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and 

(3) 

Here x Is in the direction of mean TIBL wind, y the crosswind direction, L(x*) 
the TIBL height, H(x') the plume rise, and ^^ (x 1 ) the vertical dispersion 
parameter in the stable layer. 

i 
(j (x,x') is defined as, 

where tfX. (x 1 ) is cross-wind dispersion parameter in the stable layer, and ffZ (x,x") 
the standard deviation of position of a particle at (x,0) released at (x, y', L(x')). 
Q- y Cx') and fy, (x,x/) are derived in the following sections. 

The model assumes no wind directional shear between the layered average 
wind inside the TIBL and the mean wind in the stable layer, If this shear is 
substantial then equation (1) would have to incorporate a double integral, the 
integration limits for the y' direction being a function of x'. Equation (1) is derived 
by assuming the limits on y" to be between -oeand + <p0 . The wind direction shear 
is normally between 10-20° which makes equation (1) applicable in most cases by 
assuming the mean wind direction inside the TIBL to be somewhere in between the 
direction of the plume level wind in the stable layer and the surface layer wind 
inside the TIBL. 

The rationale behind the assumption of uniform vertical distribution of 
pollutants inside the TIBL is as follows. Lamb (1978) describes the structure of the 
convective boundary layer to consist of large scale updrafts and downdrafs. The 
cores of the updrafts are more turbulent than the downdrafts. The fractional 
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area occupied by downdrafts increases with height inside the convective iayer. 
Thus particles released at the top this layer would, on the average, move towards 
the ground. Also, the magnitudes of downdraf t velocities near the top of the layer 
are small. Consequently the particles would take longer to reach ground. They 
would thus be exposed to greater number of updraf ts. This would, in turn, result in 
greater turbulent diffusion of particles released at the top. If Lamb's (1979) 
expression for within a convective boundary layer is used then the particles, 

which are diluted vigorously by the updrafts after their release at the top, would be 
expected to reach ground at a downwind distance given by, 

M 
1.65 W» y a fl x '^ 

IT 

Typically A fc£ 4 m 14 and W*/U 0.2. Therefore, equation (4) translates to a 
distance of 147 m which implies that mixing in the vertical direction is almost 
instantaneous inside the TIBL. 

A vertical distribution of the particles upto a downwind distance of about 
150m after their release at the top of the TIBL is perhaps more realistic. 
However, since an integral is used to compute downwind concentrations, the error 
resulting from neglecting a vertical distribution is not expected to be severe. This, 
together with the uncertainity in <p" inside the TIBL, which has been merely 
extrapolated from Lamb's (1979) study of a homogeneous convective boundary 
layer, does not justify the extra complexity in equation (1) which would arise from 
inclusion of a vertical distribution of pollutants inside the TIBL. Wind tunnel 
simulated thermal internal boundary layers are recommended to verify the above 
hypothesis further. 

2. RESULTS ; 

The results of the experiment are divided into two sections. First of all the 
dispersion of the plume in the capping stable layer is considered. Secondly, 
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dispersion inside the TIBL is verified with data. 

2.1 DISPERSION IN THE CAPPING STABLE LAYER; 

In general the cros-wind concentration profiles in this layer, show a near 
Gaussian shape in confirmity with model assumption. Since the data in this 
experiment were obtained at only three elevations it is difficult to assess if a 
Gaussian shape in the vertical direction is realistic. Most likely the vertical shape 
is non-Gaussian. However as discussed in section 3 assuming a Gaussian 
distribution in the vertical should not incorporate serious errors to our final 
calculation of the glc inside the TIBL. 

With this in view, the time averaged helicopter concentration data close to 
the stack were best fitted with Gaussian profiles. Following the procedure 
described in Appendix A the best values of the diffusion parameters <f*y ,(7*2 , and 
plume height H were obtained. These values are shown in Table-!. The 
corresponding theoretical estimates of these parameters are also shown in the 
same table. Here, the plume buoyancy flux Fo is obtained from the plant output 
following the procedure described in Appendix-B. 

It is seen from Table- 1, that(T y andO*Z closely follow the formulas 
suggested by Misra and McMillan (1980). Misra (198Gb)bbtained similar results for 
the data of a previous study (Portelli, 1979) at Nanticoke. The plume rise, however 
behaves in a slightly different manner. It is noted that on many occasions the 
minisonde temperature profiles show a near-neutral atmosphere at plume level. 
This may be attributed to the convective activities on the southern shore of Lake 
Erie, which can create a deep neutral layer. As this air mass is advected over the 
lake, stabilisation due to heat loss to the water surface occurs only over a small 
depth as a result of small fetch over water and moderate wind speeds. 
Consequently the stable layer capping the TIBL on the northern shore is smaller in 
depth and appears as a kink in the termperature profile. Since the level of small 



- 54 - 



scale turbulence is relatively little in the layer above the TIBL, -the plume is 
expected to rise for a long period if the temperature profile is adiabatic. Table-1 
shows the values of the Brunt-Vaisala frequency for all the data points. Employing 
Brigg's (1975) criterion for levelling off of the plume in a stable layer, it is seen 
that except for two points the plume had levelled off at the sampling location. 
Plume rise in the levelling off stage is well described by Briggs' (1975) formula. 
For the other two points where the plume was still rising, the following formula 
seems more appropriate than Briggs' (1975) formula for neutral plume rise. 

■*" * (5) 



.K - OA (£) • (JL) 



In computing plume rise only the maximum value of Fo from the two stacks is used, 
in view of the large initial separation of the stacks. The reason for the discrepancy 
in our data and Briggs formula may be expalined by the large value of entrainment 
constant which is obtained for the Nanticoke plume. It is difficult to form a 
general conclusion on plume rise at Nanticoke from only a few points. Further 
studies are required to investigate the nature of plume rise in this interesting 
situation. For the rest of the analysis we shall assume equation (5) as valid. 

2.2 DISPERSION WITHIN TIBL : 

The values of <T Yy U,x') within the TIBL are based on the dispersion 
parameters for elevated stacks in a convective boundary layer suggested by Lamb 
(1979). Thus we get, 

and (TZ^U x k) - % i*t* L*~*y 7 > (*-x') > Jb2£ 

where W* is the convective velocity scale defined as / ^Q Q U N^'HereQo is the 

a layered average temperature, and g the acceleration 
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due to gravity. Misra (1980) shows that W* is essentially constant with downwind 
distance within the TIBL. 

The shape of TIBL is assumed to be parabolic following Venkatram (1977) and 
is given as follows: 

LOO ^ Ax'* (7) 

This is well verified in a previous study at Nantiocke (Potelli, 1979). Unfortunately 
for our study a malfunctioning electronic circuitry put the acoustic sounder out of 
operation. Thus we are unable to verify (7) further. However, we believe that 
equation (7) is a reasonable characterisation of the shape of the TIBL at Nanticoke. 
This was observed also by Lui (1979). For the model the value of A is provided by 
knowing L at a given downwind distance, which is defined as the distance from the 
shoreline to the point of measurement in a direction parallel to the mean wind. In 
our case L is obtained from the inland minisonde temperature profile. 

The surface heat flux Qo was assumed to be a fraction of the net solar 
radiation. In reality Qo will vary with downwind distance inside the TIBL 
(Venkatram, 1977). However, since the convective velocity scale W* remains 
nearly constant within the TIBL, its value at a point within the TIBL sufficiently 
away from the shoreline is expected to be a representative W* for all other points 
inside the TIBL. Presumably Qo at this distant point is a function of solar 
radiation. Figure -1 shows the relationship between Qo and net solar radiation at 
such a point. These data were obtained during a previous study (Portelli, 1979). It 
is seen that the following equation represents surface heat flux reasonably well. 



Qo * 0-39 £ 



(S) 



where £\ is net solar radiation in MW.m" . We assume that Qo obtained from (3) 
is a reasonable estimate of surface heat flux at the inland minisonde site. Since W* 
is a function of (Q ) and since the model in relatively insensitive to W*/ij* 
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(Misra, 19S0 b), a slight error in Q is not expected to incorporate significant 
errors to our calculation of glc. 

Table-2 shows the parameters measured and used as input to the model for 
three fumigation events, Table-3 shows the comparison of model predicted 
concentrations to observed concentrations. It is noted that the model assumes 
uniform mixing in the vertical within the TIBL. In general the ground level 
concentrations agree within 10% with the 300' concentrations (Table-3). Thus we 
have used the 300' helicopter data to compare the model predictions on two 
occasions. It was observed in general, that helicopter and ground level Sign-x 
profiles obtained by averaging over 3-4 passes inside the plume agreed well with 
the Yz hour average concentrations obtained in a fixed mode by the mobile 
chemistry laboratory. We believe that our helicopter data at 300' of Table -3 can 
be assumed to represent the Yz hour average concentrations at corresponding fixed 
ground level monitors. 

The comparison is made on the values at different points on the cross-wind 
profiles. The cross wind profile is predicted to be relatively flat up to 500m on 
both side of its centre of mass. The observed data show a similar behaviour. This 
is to be expected as we are measuring the integrated effects of an elevated area 
source. A point source plume's cross wind profile would be much sharper. Our 
predicted profiles tend to fall off more rapidly beyond 500m from the centre of 
mass than the data. However, the concentrations levels in the tail ends being 
relatively low, are attributed a larger margin of error. Thus it is difficult to judge 
the comparison of predicted concentrations to observed concentrations at the tail 
ends. 

Figure-2 shows a scatter gram of the observed to computed concentration 
levels within the TIBL. It is seen that the data points are evenly distributed about 
the one to one line with more than 90% data remaining within a factor of two, 
suggesting that the model performs reasonably well. Similar comparison was 
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obtained by Misra (1980 a) for a previous data set. It is noted that model predicts 
ensemble average concentrtions which may not be achieved in a onveetive 
boundary layer in Yi hour. Thus the observed data based on V2 hour averaged values, 
or averaging over 3-4 passes of helicopter data, are expected to form a random set. 
It would, therefore, not be possible to get the observed values to match the 
predictions exactly. It is encouraging to note that our data remain very close to 
the one to one line, suggesting that most of our data points belong to the peak 
region of their distribution function. In air quality model validation, therefore, a 
degree of subjectivity is inevitable. 

On all occasions we observed the plume level winds to be in a different 
direction from the 85 m tower mean winds. Presumably the 85 m tower mean 
winds correspond roughly to the mean TIBL winds. Therefore, we have used the 
mean winds at the 55 m level of the tower in our model. Most of the time the 
centre of mass of the cross-wind profiles fell in a direction intermediate between 
the directions of plume level winds and the mean S5 m winds. In the model we 
have assumed that the mean TIBL wind direction is coincident with the line joining 

the source and the centre of mass of the profiles. 

2 

We observe that the plume level N plays an important role in the estimates 

of gle by affecting plume rise. Since the model is sensitive to plume rise errors in 

2 

the estimation of N at the plume level can incorporate significant errors to the 

glc. 



- 58 - 



CONCLUSIONS : 

In this paper we have shown the validation of a shoreline fumigation model 
against the data of a field study. The model is shown to reporduce the real 
situation reasonably well, Parameters describing the state of the on shore air mass 
appear to be ciritical to the final estimate of the glc. These parameters, in turn, 
depend on the convective activities on the opposite shore of the lake. This 
necessiates careful measurement of these parameters for routine application of the 
model. 

The model is simple in design and is not expected to perform well under lake 
breeze return flow situations. It Is gratifying to note that it performs as well in 
other fumigation situations. 
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APPENDIX-A 
DETERMINATION OF y . and z . 

The plume in the near field was sampled at three altitudes by a 5ignx-S0 2 
analyser installed in a helicopter. The sample air was drawn in through an 
isokinetic nozzle located on the landing skid of the helicopter clear of the down 
wash effects of the blades. The mass flow was continuously monitored during the 
flight. 

The helicopter flew an S pattern starting at the lowest level, proceeding to 
the upper levels and then descending to the lowest level. The profiles at each level 
were averaged to form an average profile to which a Gaussian distribution was best 
fitted in the following manner. 

If the concentration C inside the plume is distributed normally, then we get, 






where H is the plume's effective height Q the emission of SO~,*\/* the mean wind 
speed and Z the nth sampling height. Since, QjCy., andQTz , and H are constants 

for the entire plume during the sampling period, equation (1) implies thatlJ^C is 

2 

linearly related to y at a given sampling height. Thus if the plume is sampled at 3 

heights; viz. 2, Z~, Z?, then from the intercepts on the abscissa and ordinates of 
the best fit lines between y andx/lC at these levels we get the following: 



■ — 



where R - ■ ■ n = 1. 2, 3, and a and b are the intercepts on the 

TfVr Y «\ ' ' ' n n 

ordinate and the abscissa respectively. Equations (2) can be solved explicitly to get 
(py^Cz. and H. 
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Table 1 shows the values of y , _ , „ . . . , . ... _ „ .. „ 

J s* Z , and H obtained in this manner. Also, 

shown in this table the emission rate Q computed from the knowledge of A and 

derived from stack data. It is seen that the two Q values compare facourably. 

This shows the accuracy of the above method, 



- 62 - 



* i 



APPENDIX-B 

ESTIMATES OF STACK PARAMETERS : 

The stacks have I flues each corresponding to 4 units, each unit with a total 

load capacity of 500 MW. The gases through the flues are ejected under their 

buoyancy. Figure- 3shows the relationship between exit temperature, T , and exit 

velocity, W , of the plume v/ith load L of the power plant. It is seen that T and L 

are linearly related, whereas W versus L shows a non-linear relationship. The 

following empirical relationships are obtained from Figure-^. 

W = 11 ms" 1 , L^IOO MW 
s ' 

= 1.1 L , 100CL£ 200 MW (1) 

= 0.47L A3 , 200£I>500MW 
T s = 112 + 0.05 L, Q£Lfe 500 MW 



When n flues are operational per stack the exit velocity of the combined 
plume is obtained as 

oft * &\,+ ■-■ * dl* 

with the assumption that sufficient mixing occurs between the plumes of the flues. 
The diameter of each flue is 5.49 m. When n flues are operational the equivalent 
diameter is taken as 

d< -- \[T\ A (3) 

i.e. the areas are simply added together to from the final plume. The exit 
temperature is computed by weighting the temperature of the plume from each 
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flue by the corresponding exit velocity; i.e. 



This method assumes that total heat content of the combined plume is equal to the 
sum of the heat contents of the individual plumes. 

The buoyancy flux per stack is then computed as follows: 



Fo ^ <T" w«. (_ T fe -"*"« ■) <j 



(5) 



where Ta is the temperature of ambient air. It is seen that Fo can be expressed as 
a function of the load of the power plant. 
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TIME 


PLANT 
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. niw 


Fo 


X 
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N 
s- 1 
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iPUTED 


q 

Computed 

gs _1 

4.52 x 10 3 
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DATE 


H 

in. 


in. 


z 

in. 


II* 

m 


Cy** 

m 




from load 


May 31/79 


1400-1430 






3.7 


4.0 


380 


375 


145 


9.72 x 10 =3 




June 6/79 


1555-1630 


961 


585 


2.7 


5.0 
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262 


90 


8.26 x 10=3 
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97 


4.26 x 10? 


3.77 x 103 i 

i 


June 7/79 


1315-1340 
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2.4 


10.0 
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232 


55 





332.7 
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45 


4.67 x 10 3 


3.77 x 10 3 


June 14/79 


1400-1420 
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585 


2.7 


8,0 


393 


151 


90 


1.35 x 10^ 2 
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136 


61 


3.27 x 10 3 


3.77 x 103 

1 


June 14/79 


1655-1720 
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2.7 


8.6 


431 


219 


75 





367.5 


126 


56 


4.1 x 10 3 


3.77 x 10 3 ! 
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TABLE 2 



CASE 

NO. 



DATE 



June 13/79 



June 14/79 



June 14/79 



TIME 



1700 hrs 



1400 hrs 



1600 hrs 



ms 



-1 



7.2 



7.5 



8.61 



S -? 

wm L 



600 



1078 



814 



Qo _9 

H wm L 



240 



420 



317 



W* 



/U 



0.2 



0.25 



0.20 



-1 



1.08 x lO- 2 



1.35 x 10-2 



1.35 x 10-2 






2.4 



5.69 



5.23 



Q ys-l 



2.96 x 10 3 



3.8 x 10 3 



3.8 x 10 3 






TABLE 3 



'Oli 



CASE 


1 


2 


3 


Y 

km • 


X 
kin 


OBSERVED 
ppb (300' ) 


COMPUTED 
ppb 


X 
km 


OBSERVED 

ppb (glic) 


OBSERVED 
ppb (300* ) 

195 
195 

190 

182.5 
177.5 
167.5 
152.5 
145 
87.5 


COMPUTED 
ppb 

228 
227 

223 
217 

209 
187 
179 
148 
29 


X 
km, 


OBSERVED 

ppb (300' ) 






50 

100 

150 
200 
300 
400 

500 
1000 


17.2 


60 
62.5 

62.5 

62.5 

62.5 

50 

57 

37.5 

32.5 


48 

48 
47 
46 
45 
41 
39 
34 
TO 


16.7 


185 
185 
180 

167.5 
165 
155 
145 
135 
75 


21.3 


130 
125 
115 

100 

92:5 

82 . 5 

72.5 
35 







23-1 


110 


74.4 






— ■ 




1100 








16.77 


70 


n 





COMPUTED 
ppb 

166 

166 
163 
160 
155 
14 3 

im 
111 



Linear Regression 
r 2 = 0.85 



c 0BS = 28.83 +0.65 C C0Mp 
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Appendix II 
Daily Flight Description During Nanticoke May/June 1979 Field Study 
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MAY 30 DATA 
FLIGHTS 5 AND 6 
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May 30 



May 30 was a generally sunny day with some cloudy periods. The winds 
were light from NNW in the morning, shifting to SSW by late afternoon. The 
ground based units began to report ground level SO- at 1415 hrs in the HWY 6 
and HYW 3 area. The helicopter began sampling at 1445 hrs. 

Flight 5 

Flight 5 began at 1445 hrs, at which time the plume was found to be aloft 
over Jarvis at 1200 feet above ground level. The first ground level SCU was 
detected near Hagersville, and a set of traverses were attempted at this point. 
However, part way through the first set of traverses, the plume became very 
dilute, and the set could not be completed. Another set of traverses was 
attempted upwind of Hagersville, howerver, no significant ground level SO- was 
found. The flight was completed with a set of traverses near the stack. 

The data collected in flight five were of very little use because a complete 
set of traverses was not completed in the fumigation area, and as such the data 
is not presented. 

Flight 6 

After refueling the helicopter began to sample for flight 6 at 1735 hrs in 
the Hagersville area along HWY 40. Six traverses were completed ( 2 at 300, 600 
and 900 ft) in this area, during which time the plume was found to be uniform 
with height. However, the levels detected were quite low, and by the time the 
sixth pass was completed there was very little SO- detected. Flight six was 
completed with a set of traverses near the stack, along Concession 1. Six passes 
were completed, 2 at each altitude of 1200, 1400 and 1600 feet. 

Filter packs were exposed near the stack at 1200, 1400 and 1600 feet above 
ground and in the fumigation zone at 300, 600 and 900 feet. Hydrocarbon 
cartridges were exposed in the fumigation zone, and up-wind, over the Port 
Dover area, 
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MAY 31 DATA 



FLIGHT 7 
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May 31 



May 31 was a sunny day, with some cloudy periods. The winds were from 
east throughout the day. The COSPEC unit indicated that the plume was running 
parallel to the shore most of the day, however, it was over Port Dover and 
fumigating at 1230. 

Flight 7 

The helicopter began sampling for flight 7 at 1255 hrs NW of Port Dover, 
however, upon arrival the plume had shifted over the lake. The plume was then 
followed to the Turkey Point area where it was found to be fumigating. At this 
time the plume appeared to be very brown towards the stack. Two passes were 
completed at 300 and 600 feet, while only one could be completed at 900 and 
1200 feet before the plume started to move out over the water and become more 
coherent. 

The plume was then followed back to the chimney for a set of traverses 
over the STELCO pier. Eight traverses were completed in this area, with 2 at 
900, 1100, 1300 and 1500 feet. 

The day's sampling was completed with a background flight over Peacaock 
Point at 1100 feet. 

Filter packs were exposed at all three locations. The data collected is of 
marginal use because it was not possible to repeat the traverses in the 
fumigation area. 
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FLT 7/2 
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FLT 7/2 



1434 HRS 
1500. FT 
PASS 7209. 
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JUNE * DATA 
FLIGHTS 8 AND 9 
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June 4 



On this day units i, 6 and 7 at the generating station were operational and 
the weather was sunny with southwest flows. 

FLT 8 The ground based units spotted fumigation at 1330 and the helicopter 

began sampling at 1357 along Hwy 3 between Nells Corners and Cayuga. The 
plume was visible upwind, but fumigation was taking place at this position as 
seen by the high concentration at 300 ft a.g.l. A traverse was carried out at a 
bearing of 90° and as such may not have been crosswind, however, dust blowing 
from a local quarry suggested this was the best direction. The aircraft speed was 
quite variable around 30 mph (20-40) with the average speed around 30-35 mph. 
The plume was patchy and difficult to locate, and after passes at 300, 600 and 
900 ft. the plume shifted and could not be relocated. The flight was completed 
by making 3 passes 1 km NE of the stack where the plume was found to be at 
about 1200 ft a.g.1. Filter packs were exposed in the fumigation area and at 
stack, and a hydrocarbon sample was taken in the fumigation area. 

This flight seems to provide marginal data in that only three passes were 
made at each down-wind distance, however the MAMA and COSPEC units were 
sampling at the same time in the same area. There is also some question 
regarding the bearings chosen for the flight traverse because the wind directions 
from pibal data show the wind direction to be approximately 240 from N. 

FLT 9 After refuelling, the helicopter began sampling for flight 9 at 1635. 

The mobile units were sampling in the Hwy 40 area at Nells Corners, however, 
the plume was very visible further downwind. The first significant SO- levels at 
300 ft a.g.l, were found just east of the Grand River and a flight path at bearing 
150° was established between endpoints H2 and 12. The wind was gusty and 
therefore it was difficult to maintain constant bearing and speed. 

While sampling in this area the plume was light brown and visible on all 
sides. A very dark area was visible further downwind. 1 The plume in the area 
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was well defined but seemed to be in layers. Six passes were made in the area (2 
at 300, 600, adn 900 ft) and much of the plume seemed to be aloft. After six 
passes were completed the helicopter moved about k miles downwind to sample 
in one of the concentrated areas of the plume. This sampling was between 
endpoints L2 and K2 at a bearing of 290°. Three passes were made at this 
location (at 300, 600 and 900 ft) and there still seemed to be more plume aloft 
than at the lower levels. 

The plume was then followed in to the stack and traverses were made at 
1200, 1000, 800 and 600 ft with the plume height at approximately 1100 ft. 

Filter packs were exposed in all three areas and the flows monitored. The 
Sign-X baseline was stable throughout the day. 
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JUNE 6 DATA 
FLIGHT 10 
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June 6 



June 6 was a mainly sunny day with some hazy periods. The winds were 
generally from the east in the morning, shifting to east -south -east early in the 
afternoon. 

No Sign- X calibration was performed, therefore an average value for the 
slope of the calibration curve was used. The baseline and flow were stable all 
day. 

Flight 10 

The ground based units began to report ground level SO ? at 1430 hrs north 
of Port Dover along Road 5. The helicopter began sampling in this area along a 
flight path between endpoints S2 and X2 at a bearing of 30°. The aircraft speed 
varied between 20 and 40 mph, averaging 30 mph. 

A total of 10 traverses were completed at altitudes ranging from 300 to 
1300 feet above ground level during which the plume was found to be spread out 
evenly along the entire length of the traverse, and homogeneous with height. 
During the sampling period the plume was shifting rapidly towards the north as 
the lake breeze component of the wind became stronger. 

The flight was completed by flying eight passes at altitudes of 2200, 2000, 
1800 and 1600 feet over the Nanticoke Raod - Concession 1 area (2 km from the 
stack). 

The day's sampling was completed by flying a background flight (Fit 1 1) for 
30 minutes over Peacock Point. 

Filter packs were exposed in all three areas, and hydrocarbon cartridges 
were exposed in the fumigation zone and during the background flight. 

The data is very useful not only because of the large number of profiles 
collected in the fumigation area, but because the COSPEC vehicle was sampling 
in the same area during the sampling period. 
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JUNE 7 DATA 
FLIGHT 12 
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June 7 



June 7 was a warm, mainly sunny day with winds from the south-southeast. 
It was hot and hazy in the morning, with cloud developing in the afternoon. 
During the helicopter sampling period it was generally hot and humid with 
variable cloudiness. The ground based units began to report gound level SO- at 
1300 hrs, along Highway 40. 

The helicopter SIGN-X was calibrated before the flight, however, periodic 
flow adjustments were required during the flight, causing changes in the 
baseline. 

Flight 12 

The helicopter took off at 1300 hrs and began traversing near the stack in 
the Sandusk Road - Concession 1 area. Six traverses were completed in this 
region ( 2 each at 800, 1000 and 1200 ft) which determined the plume height to 
be 1000 feet. 

The helicopter followed the plume downwind and located the fumigation 
area near Highway 40. Sampling was then carried out paralled to Highway 40 
between Sandusk Road and Fisherville Road. Six traverses (2 at each altitude of 
300, 600, 900 ft) were completed in this area before the fumigation area moved 
further downwind. The plume was well defined, and visible during this run. 

The helicopter then moved further downwind and found the fumigation area 
to be near the Grand River (approximately 10 km NE of Highway 40). A set of 
traverses was then completed in this area between endpoints 13 and E3 (6 passes, 
2 each at 300, 600 and 900 ft). 

The sampling was completed by making a background flight in the Port 
Dover area (Flight 13). 

Filter packs were exposed during the near stack passes, in each of the 
fumigation areas, and in the background flight. Hydrocarbon cartridges were 
exposed in run 2 (along Highway 40) and during the background flight. 
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June 13 

June 13 was a warm, sunny day with some cloudy periods. The winds were 
steady, from the west -south west throughout the day. 

The ground based units were sampling along Highway 40, and began to 
report ground level SO- at 1500 hrs, The helicopter took off at 1515 to begin 
sampling for Flight 18. The SIGN-X was calibrated prior to the flight, and the 
baseline and flow remained constant all day. 

Flight 18 

Sampling began in the suggested area along Highway 40 between Highway 3 
and Rainham Road, however, after two passes there were no reasonable levels so 
sampling was attempted in the Rainham Road, Highway 3 area - again with no 
success. The helicopter then moved closer to the stack to find the plume, and 
then followed it downwind to an area NE of Highway 40 between Nells Corners 
and Fisherville Road. Samping was then carried out between endpoints 34 and 14. 
The plume was found to be quite "patchy" and spread out evenly over the length 
of the traverse. Three passes were completed at each of the three altitudes 
(300, 600 and 900 ft). The flight was completed by making a set of passes near 
the stack in the Concession 1 - Selkirk Road area. Twelve traverses were made, 
with three at each altitude (600, 1000 and 1200 ft). The plume height was 1100 
ft. 

The helicopter sampling for June 13 was completed with a background 
flight (Flight 19) up-wind of the stack in the Port Dover area at 1100 feet 
above ground level. 
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June 14 

This was a warm sunny day with south -south west flows which were light in 
the morning increasing in the afternoon. Two units were operational at the 
generating station. Two flights were performed this day, and all ground based 
units were operational. The Sign-X was calibrated in the morning and the 
baseline remained reasonably constant throughout the day. 

FLT 20 

The mobile monitors started to pick up significant ground level 
concentrations at about 1242, and the helicopter was in the air sampling for 
flight 20 at 1300 hours. The flight path for this flight consisted of a traverse 
paralled to Hwy 40 between endpoints B5 and F5 (6.5 km with bearing of 110°). 
A total of 15 passes were made at three altitudes (5 passes each at 300, 600 and 
900 ft). The wind seemed to be veering during the sampling period as the plume 
was moving eastward. The plume was distributed evenly in the vertical direction 
as seen by the constant levels measured at 300, 600 and 900 ft. and was quite 
brown and visible for a long distance down-wind. The aircarft speed was varying 
between 20-40 mph and the aircraft was being blown down-wind which 
necessitated frequent adjustments in the bearing of the aircraft. 

A set of passes were also carried out at the stack. Nine passes were made 
at 3 altitudes (1000, 1200, and 1400 ft) and the plume height was found to be 
approximately 1200 ft. This flight was completed with a background flight west 
of Nanticoke. 

Filter packs were exposed at all three locations and hydrocarbon cartridges 
in the fumigation and upwind areas. 

This flight provides very good data because of the large number of passes 
(5) at each altitude in the fumigation area and near the stack. The helicopter 
was also sampling in the same area as the ground based units which will provide 
good data for comparison. This was a very good example of fumigation. 
FLT 21 

Sampling began at 1545 between endpoints N5 and Q5 at a bearing 
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of 115 located about 5 km NE of Hwy 40. Wind was reasonably stable during this 
run and as such the plume did not shift much. The helicopter speed varied 
between 20-40 mph, with the average speed being 30-35 mph. The plume was 
again found to be evenly distributed vertically, and the fumigation area was 
clearly visible from the helicopter. Twelve passes were made at this location 
with 4 passes at each of 3 altitudes (300, 600 and 900 ft). Part way through the 
sampling there was some blasting occuring at a gravel pit upwind of the 
helicopter, and therefore some particulates may have been sampled. There was 
not noticable effect on the Sign-X trace, and the sample line was clean when 
checked after the flight. 

Eight passes were made at the stack with 2 passes at each of the four 
levels (1000,1200,1400 and 1600 ft)and the plume was determined to be at 1200 

ft. a.g.l. 

Filter packs were exposed up and downwind, and two hydorcarbon tubes 
were exposed in the fumigation area. 

Good data has been obtained from the helicopter for this flight due to the 
large number of traverses in the fumigation area. However, the helicopter 
sampling is located 5-6 km from both of the ground based units. The angles of 
the traverses are also quite different for the COSPEC vehicle relative to the 
helicopter. 
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Appendix III Minisonde Measurements 
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This Appendix consists of profiles of temperature, wind direction and wind 
speed as functions of height obtained from minisonde flights at two locations. 
The hourly profiles are presented in pairs, with the MSI site located at the 
lakeshore, and the MS-2 site at some distance inland, The position of the MS-2 
site varied depending on the wind direction, The positions of the various release 
sites are indicated on a map for each day contained in Appendix II. 
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